Electronic Structure and Quasiparticle Band Gap of Silicene Structures 
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We report first-principles results on the electronic structure of various silicene structures. For pla- 
nar and simply buckled silicenes, we confirm their zero-gap nature and show a significant renormal- 
ization of their Fermi velocity by including many-electron effects. However, the other two recently 
proposed silicene structures exhibit a finite band gap, indicating that they are gapped semiconduc- 
tors instead of previously expected Dirac-fermion semimetals. Moreover, our calculated quasiparticle 
gap quantitatively explains the recent angle-resolved photoemission spectroscopy measurements. In 
particular, the band gap of the latter two structures can be tuned in a wide range by applying strain, 
giving hope to bipolar-devices applications. 



Graphene, a layer of hexagonal carbon lattices, has 
spurred tremendous interest because of its unique lin- 
ear energy-momentum dispersion and exciting applica- 
tions. p]-[5| Its huge success has also motivated signif- 
icant efforts to look for the similar honeycomb struc- 
tures but made of other group IV elements, such as sil- 
icon and germanium, which are so called silicene and 
germanene, [rjl-[l7| respectively. Particularly because of 
the compatibility with the matured silicon technologies 
and the stronger spin-orbit coupling for potential topo- 
logical insulator candidates, [l3| silicene is expected to be 
a promising material both theoretically and practically. 

Because of the strong trend of electrons of silicon 
atoms to form the tetrahedral sp 3 hybridization, unlike 
graphene, a simply buckled silicene structure is predicted 
to be more stable than the perfectly planar one. [8] This 
buckling provides a unique degree of freedom to silicene 
atomistic structures; to date several unique structures 
with more complicated buckling styles are proposed by 
recent experiments, in which a layer of silicene has been 
successfully fabricated on metallic substrates. p^34l7| The 
Dirac-fermion energy-momentum dis per sion is observed 
in these silicene structures as well. jl4l-[l8j However, other 
than interpreting atomistic structures, very limited at- 
tempts have been carried out to discover the electronic 
structures of these recently proposed silicene structures, 
which is crucial for guiding experimental measurements 
and understanding the electric and optical properties of 
these materials of ever-growing interest. 

Moreover, the buckling of silicon atoms may break the 
symmetry of the honeycomb lattice. According to the 
history of graphene, the variation of structure and sym- 
metry shall not only change the Fermi velocity but also 
possibly generate a finite band gap, promising a crucial 
advantage over graphene for broader applications, such 
as bipolar devices and high-performance field-effect tran- 
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sistors (FETs). 

In this Letter, we present our first-principles calcula- 
tions, using both density functional theory (DFT) and 
many-body perturbation theory (MBPT), to reveal the 
electronic structures of several promising silicene struc- 
tures. Our study shows that both the planar and simply 
buckled silicene structures have a gapless linear energy- 
momentum dispersion and that the self-energy correction 
contributes to a significant enhancement of the Fermi ve- 
locity (~ 37%). However, for the two recently proposed 
silicene structures, a finite band gap is identified and this 
gap value can even be varied in a wide range by tuning 
the strain. The inclusion of many-electron effects by the 
GW approximation further confirms this finite band gap, 
which corrects the previous key interpretation of these sil- 
icene structures as massless Dirac-fermion semimetal. At 
the same time, our GW calculated band gap satisfactorily 
explains the recent angle-resolved photoemission spec- 
troscopy (ARPES) measurements. On the other hand, 
the tunable band gap of silicene according to the applied 
strain predicts promising device applications, overcoming 
the known trouble of gapless graphene. 

For our DFT simulations, both local density approx- 
imation (LDA) and generalized-gradient approximation 
(GGA) are applied to make sure our conclusions are not 
sensitive to the choice of functionals. [l9|, [20| The Kohn- 
Sham equation is solved using the plane-wave basis with 
a 24 Ru energy cutoff. Norm-conserving pseudopoten- 
tials [2l| are applied and the k-point sampling grid is set 
to be dense enough to obtain converged DFT eigenvalues. 
In order to provide more quantitative results with many- 
electron effects included, the single-shot Go Wo approxi- 
mation is employed to calculate QP band gaps [22|, [23| 
with a layered Coulomb truncation. The k-grid sampling 
of such GW calculations is 64 x 64 x 1 for planar and 
simply buckled silicene structures while 32 x 32 x 1 and 
16 x 16 x 1 for those proposed by recent experiments, 
respectively. 

Here we consider four typical silicene structures of cur- 
rent fabrication interest, which are shown in Fig. [TJ The 
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FIG. 1: (Color online) The ball-stick model of studied sil- 
icene structures, (a) the planar silicene; (b) the simply buck- 
led silicene; (c) buckled silicene with a (V3 x v^) unit cell; 
(d) buckled silicene with a (4 x 4) unit cell. The buckling 
distance is labeled by d. For each structure, we present both 
the top view and side view. The unit cell is illustrated by 
the shadow region. The different colors of silicon atoms are 
assigned according to their vertical positions. 



first two are known planar and simply buckled struc- 
tures while the last two are proposed by recent experi- 
ments with (4x4) and (V3x y/3) structures, respectively. 

main differences between these structures are 
the size of the unit cell and the way to buckle the lattices. 
For example, the silicon atoms in the (\/3x y/S) structure 
shown in Fig. [T] (c) are located in three horizontal planes 
while those in the (4 x 4) structure shown in Fig.[T](d) are 
located in two horizontal planes. According to the total 
energy of DFT, the simply buckled structure shown in 
Fig. [D (b) is the most stable configuration with the in- 
plane lattice constant of 3.81 A and a buckling distance 
d =0.40 A, which are in agreement with previous results. 

It has to be pointed out that those structures shown in 
Figs.[T](c) and (d) do not exactly mimic the experimental 
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FIG. 2: (Color online) (a) DFT and GW calculated band 
structure around the Dirac cone of simply buckled silicene. 
(b) The quasiparticle energy vs. LDA band energy of sim- 
ply buckled silicene, which shows the linear correction with a 
slope of 1.37. 



cases, in which the substrate may be an essential ingre- 
dient [3, This is indicated by our DFT calculations, 
and other published works show that these isolated struc- 
tures are only met ast able. [16| Therefore, as what have 
been done before, [29] our solution is to fix the buckling 
distance d and relax all other degrees of freedom, includ- 
ing force and stress. Our following calculations will show 
that the essential physical picture and main conclusions 
shall not be affected by this approximation. 

For the planar and simply buckled silicene structures 
shown in Figs. [1] (a) and (b), both DFT and GW calcula- 
tions confirm the zero-gap nature with a linear massless 
Dirac-fermion dispersion. Since planar silicene has nearly 
the same energy-momentum dispersion (with slightly dif- 
ferent Fermi velocities) and GW-LDA correction as the 
simply buckled silicene, we only provide the results of 
the simply buckled case in Fig. Fig. [2j Figure [2] (a) 
shows the band structures around the Dirac Cone for the 
simply buckled silicene, and Fig. [2] (b) shows its quasi- 
particle energy vs. LDA band energy relation. A perfect 
linear relation is observed, which means the zero-gap lin- 
ear band dispersion is kept even including many-electron 
effects. From these linear energy- momentum dispersions, 
we get the corresponding Fermi velocity of the planar 
silicene as vp = 5.6xl0 5 m/s and that of simple buck- 
led structure as Ui? = 5.4 xlO 5 m/s at the DFT level. 
After including the electronic self-energy corrections by 
the GW approximation, these values are enhanced to be 
7.7xl0 5 m/s and 7.4 xlO 5 m/s, respectively. These Fermi 
velocities are substantially smaller than that of graphene 
(~ 1.1 xlO 6 m/s) but the renormalization of the Fermi 
velocity by self-energy corrections is similar to that of 
graphene, which is around a 37% enhancement from their 
DFT result [HI, because of the depressed screening ef- 
fect. 

Then we turn to the recently proposed structures from 
experiments. In Fig. [3] (a), the DFT-calculated band 
structure of the buckled (\/3 x \/3) structure is obtained 
by using the buckling distance d = 0.71 A. which is a 
typical buckling value. [H, [H], EH, 03 Surprisingly, there 
appears to be a finite band gap (~ 0.21 eV) at the V point, 
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FIG. 3: (Color online) (a) Band structure of the (y/3 x y/3) 
buckled silicene structure shown in Fig. 1 (c) with the buck- 
ling distance as d — 0.71 A. (b) Band structure of the (4 x 4) 
buckled silicene structure shown in Fig. 1 (d) with the buck- 
ling distance as d — 0.75 A. The top of valence band is always 
set to be zero. 



which is qualitatively different from previously claimed 
massless Dirac-fermion dispersion. pj| Actually, we can 
obtain a well-defined effective mass of the free carriers 
from the curvature of the band dispersion in Fig. [3] (a). 
For example, the nearly isotropic effective mass of the 
light electron is 0.083 mo and that of the hole is 0.079 
mo, which are typical values of semiconductors. 

Moreover, we further tune the buckling distance d 
while maintaining the supercell geometry to obtain the 
evolution of the band gap, which is presented in Fig. [4] 
(a). From another point of view, this is equivalent to 
tune the strain condition of the layer structure. For both 
LDA and GGA results, the band gap of such a (y/3 x y/S) 
structure increases while enlarging the buckling distance. 
In particular, as the buckling distance is more than 0.8 
A, we see a transition of the band gap from the direct 
one to the indirect one, which is reasonable because the 
band structure shall approach that of tetrahedral silicon 
due to the stronger sp 3 hybridization. Therefore, unlike 
the planar and simply buckled silicene, whose band gaps 
are always zero under uniaxial or biaxial strain, the band 
gap of this (V^ x a/3) silicene structure can be tuned from 
zero to a wide range by the mechanical strain, which shall 
be of broad interests. 

To further confirm this finite band-gap nature, we have 
performed the GW calculation. For the purpose of justifi- 
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FIG. 4: (Color online) The band gap evolution as the vari- 
ation of buckling distance d for (a) the (y/3 x y/3) silicene 
structure and (b) the (4 x 4) silicene structure. The direct 
band gap values are connected by dash lines while the indi- 
rect band gap values are connected by solid lines. The GW 
results are marked by cross signs. 



cation, one GW calculation of a typical buckling distance 
shall be enough; as shown in Fig. [4] (a), the QP band gap 
is around 0.5 eV for the d = 0.71 A buckled case, which 
is almost a 140% enhancement from the DFT result be- 
cause of the substantially depressed screening. It must 
be addressed again that our referenced experiments use a 
metallic substrate, which may reduce the self-energy cor- 
rection due to the metallic screening and possible charge 
transfer. However, all of these factors will not close the 
finite band gap and, usually, the realistic QP band gap 
shall still be slightly larger than the DFT result. [HI [23 

We now turn to the band structure of the other impor- 
tant silicene structure (4x4), as shown in Fig. El (b). For 
this superstructure, the buckling distance d = 0.75 A is 
chosen from the experimental study. [14] This time we 
observe a finite band gap again, which is a typical direct- 
band-gap semiconductor as shown in Fig. [3] (b). The 
evolution of the band gap with the buckling distance d 
is also presented in Fig. 2] (b). Similar to the (\/3 x y/3) 
case, the band gap increases as we enlarge the buckling 
distance. When the buckling distance is more than d = 
0.8 A, the similar direct-to-indirect band gap transition 
is observed in Fig. 2] (b). The corresponding GW results 
are marked as well, predicting a larger QP band gap than 
the DFT results. 

Interestingly, when we consult the relevant experimen- 
tal result, a finite band gap was observed, but it had been 
attributed to extrinsic factors, such as the layer-substrate 
interaction. [14[ However, our simulation shows that this 
finite band gap may be intrinsic if the silicene sample 
possesses the claimed (4 x 4) structure. Moreover, if we 
compare the experimentally observed band gap value by 
ARPES to our first-principles GW result shown in Fig. [4] 
(b), they are very close to each other, 0.56 eV from our 



4 



(a) (b) 




FIG. 5: (Color online) Top view of color contour plots of the 
charge distribution of typical conduction and valence states 
at energy extrema around the band gap for (a) the planar 
silicene when d = and (b) the (4 x 4) silicene structure 
when d = 0.75 A. 

calculation while 0.6 eV from experimental work, for the 
(4 x 4) structure with a buckling distance of d = 0.75 A. 

After presenting first-principles results, it is necessary 
to figure out the physical reason for such a finite band gap 
opening in these promising silicene structures. According 
to past intensive studies on graphene, a necessary condi- 
tion for the massless Dirac cone with a zero-gap character 
is the inversion symmetry of AB sublattices. However, 
for the structures shown in Figs. Q] (c) and (d), A and 
B sublattices are no longer equivalent to each other, re- 
sulting in the broken inversion symmetry and, thereafter, 
generating a finite band gap. A similar idea was noticed 
before, e.g., the applied gate electric field can break this 
symmetry and introduced a finite band gap in simply 
buckled silicene and germanene. [28l430j . 

In order to better see this broken AB sublattice sym- 
metry, we present the charge distributions of valence and 
conduction states at energy extrema around the band 
gap. As shown in Fig. [5] (a), for planar silicene, we ob- 
serve the charge distributions of both valence and con- 
duction states are nearly identical for A and B sublat- 
tices. However, for the (4 x 4) silicene structure, we can 
clearly see the inequivalent charge distribution of A and 
B sublattices as shown in Fig. [5] (b) due to the broken 
symmetry. 
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All the above calculations consider only free-standing 
silicene structures. However, the predicted finite band- 
gaps in the structures, as shown in Figs. Q] (c) and (d), 
will persist even if the substrate effect is included be- 
cause the layer-substrate interactions shall further break 
the symmetry and enlarge the band gap. [31( In this 
sense, we believe that the recently proposed silicene 
structures have qualitatively different electronic struc- 
tures and band gaps from the widely assumed massless 
Dirac-fermion dispersion. The most likely reason for this 
conflict between our simulation results and previous ex- 
perimental conclusions may be due to the interpretations 
of atomic structures from their experimental data. 

Finally we have to point out that the finite-gap na- 
ture of these recently proposed silicene structures shall 
be of practical interest. As shown in Figs. [3] and |4j the 
band gap and even the band topology can be tuned in 
a wide range by the buckling distance d and the associ- 
ated strain; this provides a precious degree of freedom to 
control the electronic and optical properties of such ma- 
terials, making them superior to gapless graphene. Ac- 
cording to previous studies, H, GH 0, [H, LL7| the lattice 
structure and buckling distance are able to be varied by 
layer- substrate interactions. Therefore, if future experi- 
ments can fabricate and confirm more similar structures, 
they shall be of broad interest for devices applications 
with such a tunable band gap. 

In conclusion, we employ first-principles simulations to 
study the electronic structure of various silicenes. Our re- 
sults show that the band structure of silicene is sensitive 
to the buckled atomistic structure. Unlike the planar and 
simply buckled ones, these silicene structures proposed 
by recent experiments exhibit a finite band gap, making 
them promising 2D semiconductors instead of zero-gap 
semimetals, as were previously assumed. On the other 
hand, this finite band gap can be tuned in a wide range 
by strains and makes these silicene structures superior 
to graphene for many important applications. Our pre- 
dictions reinterpret relevant experimental measurements 
and shall motivate more reliable justifications. 
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